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Abstract 

Heavy multi-quark potential in the SU(N) color group using hard-wall AdS/QCD at both zero 
and finite temperature is studied. A Cornell-like potential is obtained for baryons and other exotic 
configurations and compared with those in the quenched lattice calculation in N = 3 case. At the 
end we also discuss possible improvements in the UV region of potential. 
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I. INTRODUCTION 



When the mass of a heavy quark rriQ is much larger than the QCD scale Aqcd, the heavy 
QQ bound states (quarkonium) are like hydrogen atoms where ordinary quantum mechanics 
can be applied. In the non-relativistic limit, it is possible to describe the interaction between 
Q and Q in terms of a local potential E(L), where L is inter-quark separation. By solving 
the Schrodinger equation in three-dimensions, the quarkonium spectrum can be adequately 
characterized by E(L) and itlq. This potential approach to quarkonia has been applied to 
the study of Charmonium and Bottomonium states and a fair agreement with measurement 
was found. One of the most popular models for E(L) is the so-called Cornell potential[l|, 

E(r) = -j + aL. (1) 

The first part is Coulombic due to the one-gluon exchange between Q and Q. The second part 
is the confinement potential due to the formation of color flux tube. Relativistic, radiative 
and other non-perturbative corrections can also be added ad hoc as extra terms. Though 
we are still lacking a first-principle derivation of the whole potential, among few approaches 
to the non-perturbative QCD, lattice QCD computation provides an irreplaceable way to 
test its validity. 

On the other hand, the color flux tube can be easily modelled by a string. In particular, 
the coefficient a in the equation ([1]) is identified as the string tension. Though this old 
version of String Theory first failed to be a realistic QCD model for its false prediction of 
some massless particles, it nevertheless took another ambitious task as a candidate of Theory 
of Everything and later evolved into its supersymmetric counterpart, the Superstring Theory. 

However, not only QCD string survived but in fact it revived by gaining new insights 
from the so-called Anti-de Sitter space-Conformal Field Theory (AdS/CFT) correspondence, 
first proposed in 

Baa. 

In particular, it was applied to relate the thermodynamics 
of M = 4 super Yang-Mills (SYM) theory in four dimensions to the thermodynamics of 
Schwarzschild black holes in five-dimensional Anti-de Sitter spacej^]. In this description, 



confinement /deconfinement phase transition of gauge theory on a sphere has its holographic 
dual description as the Hawking-Page phase transition p]. Later, it was realized that con- 
finement can be achieved by capping off the Calabi-Yau cone smoothly at the infrared 
tipQ, S 3, or by introducing IR cutoff in the AdS space 3, 11, 12, 3, 14, 15, Q, Q. 



This latter approach are usually referred to as the bottom-up construction of AdS/QCD, 



together with the other top-down construction by [la, [19(, are two main approaches to 
realize QCD physics via the (super)gravity theory in specific backgrounds. In particular, 
it is hoped that, as an alternative to lattice QCD, the non-perturbative region of QCD 
can be better understood both qualitatively and quantitatively within its dual picture of 
gravity. One important testing ground for the AdS/QCD is the experiments of Relativistic 
Heavy Ion Collision (RHIC), collisions of gold nuclei at 200 GeV per nucleon are about to 
produce a strongly-coupled quark-gluon plasma (QGP), which behaves like a nearly ideal 
fluid. While the perturbative calculation can not be fully trusted in this strongly-coupled 
region, there are increasing amount of interests in calculation of hydrodynamical transport 
quantities via the use of AdS/CFT correspondence, in particular that the energy loss of a 
heavy quark moving through M = 4 SYM thermal plasma has been extensively studied 20] . 
The drag force was derived in the context of AdS / CFT to model the effective viscous in- 



teraction 
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, l22j, later it was generalized to a rotating black hole or with a dilaton field 



26j and B-field 



27 



281 ]. Drag force of a comoving (QQ) pair was also considered 



291 ] and energy loss of baryon was studied in 30] . It is hoped that this line of research 
will eventually make contact with experimental results from RHIC. 

Back to the equation (JT]), one would like to know how to derive it in the revived String 
Theory. We have learnt that (QQ) static potential can be calculated via temporal Wilson 
loop for time T>i, 



< W >~ e 



-TE(L) 



(2) 



In the AdS / QCD scenario, this quantity, at zero genus, can be calculated via on-shell classical 
(super)gravity action in the AdS bulk geometry, i.e. the minimal surface enclosed by the 



Wilson loop on the boundary 



311 ] . Similarly, baryon potential can also be constructed once 



the baryon vertex is realized as a wrapped brane on the compactified sphere 32j, 



34 



35, 



361 ] . However, the potential obtained in the AdS background is always Coulombic thanks 



to its conformal symmetry. A Cornell-like potential for heavy (QQ) was obtained in Jjj] by 
breaking this symmetry via new scale set by an IR cut-off. This so-called hard-wall model 
simply puts an IR cut-off brane at location r = R. This cut-off brane is responsible for 
confinement as it will become clear later. The corresponding metric is then, 



ds 2 = (-) 2 (-dt 2 + dx 2 ) + (-fdr 2 + R 2 dQ 
R r 



(3) 



where R < r < oo. 
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In this paper, we go beyond the study of heavy (QQ) static potential, and look for 
multi-quark potential in the SU(N) color group using the above-mentioned toy model of 
AdS/QCD at both zero and finite temperature. 

This paper is organized as follows. In section [Til we refresh the construction of heavy 
(QQ) potential as in [ljj]. In section UTTI we construct the heavy baryon potential and make 
a naive comparison with that in the quenched lattice calculation for SU(3) color group . 
In section IIVI we construct exotic multi-quark configuration. In particular, we obtain tetra- 
quark and penta-quark potential and compare them with the lattice results. In section [V] 
we study baryon potential at finite temperature. In section |VTJ we make some proposals to 
improve our construction, but also discuss their limitations. In section IVIIt we conclude 
with summary and a few comments. 



II. HEAVY MESON POTENTIAL 

In this section, we recall the construction about inter-quark potential of heavy (QQ) in 
truncated AdS space. When the quark and anti-qaurk are close enough, the potential E 
and inter-quark distance L are given by 

*M = ^(f *<^=-i)-i), 

p2 poo i 

L(r )=2- dy (4) 
r J i y 2 ^y 4 - 1 

where r is the lowest point where QCD string can reach in the bulk. The AdS (S 5 ) radius 



R A = A/a' 2 = gy M N/a' 2 . Therefore one obtains Coulomb-like potential 31] 



E = a ~ 0.228. (5) 

This result is expected from conformal invariance. As the distance increases, the U-shape 
string will be cut by the IR brane and the potential will include an additional term, given 



by 



E' = E(R) + ^- ; (L-L(R)). (6) 



As a result, we recover linear potential from the second term for large separation. 
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FIG. 1: Inter-quark potential as a function of baryon size. Fitting curve A agrees with the lattice 
simulation at IR region, while curve B agrees at UV region. 

III. HEAVY BARYON POTENTIAL 

The baryon (Q ■ ■ - Q) in the SU(N) SYM is described by iV open strings joining at a 
vertex formed by wrapping a D5 brane on S 5 . The introduction of open strings is essential 
to incorporate quarks in fundamental representation of the SYM. The total energy of this 
baryon configuration is simply a sum of energy of N Fl's and a wrapped D5 brane, given 



where the vertex is situated at the center of AdS§ and the total energy has been regularized 
by subtracting off the energy of free quarks. 

With the IR cut-off, we simply allow the configuration of each open string to be the same 
as in the original AdS space as if the vertex were at r = ro < R. However, those strings 
are cut by the IR brane and replaced by string segments on the brane. Therefore the vertex 
is actually projected onto the brane. The relation between r and baryon radius L is still 
valid, i.e. 



by 33] 




L 



(3 ~ 0.007 



(7) 




7 



0.481 



(8) 



and the projected radius Lr on the cut-off is 




i yVW-l' 



dy 



(9) 



where y = y R = ^ and (3 = a/16/15 



The total energy, composed of one vertex and 
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FIG. 2: Flux-tube recombination between connected four-quark state and two-meson state is re- 
ferred to the flip-flop. It happens as two vertices are close enough. 



strings, becomes 




(10) 



After replacing r and E in terms of L, and choosing properly, it becomes 



E = —N— + aNL + C. 



(11) 



In particular, at the limit tq 




i 



A 



1 



C 



R 5 



(12) 



2na' 




In the FIG. 1, we compare with the quenched lattice result for 577(3) color group[37|, i.e. 
(A,a,C) = (0.0768,0.1524,0.9182). However, two variables are generally insufficient to fit 
three unknowns. If we fit the same a and C with the choice of a' = 1.0443 and R 2 = 2.2591 
(or A = 5.1036), then a twice bigger A = 0.1656 is obtained for curve A. If we fit the same A 
and a, then we obtain C = 0.1345 (or A = 1.0064) for curve B. As a result, fitting curve A 
agrees with the lattice simulation in the IR region, while curve B agrees in the UV region. 

IV. EXOTIC MULTI-QUARK POTENTIAL 

One can imagine that multi-quark system might be formed in the process of pair anni- 
hilation as two or more baryons/anti-baryons are dragged toward each other. In the case 

1 L here refers to the projected distance from constituent quark to the vertex, not the inter-distance between 
quarks as used in the lattice simulation [3?]]. If the convention of lattice is used, we obtain A' ~ 0.133\/A, 
about half the coefficient of meson. 
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FIG. 3: To the left: It shows that penta-quark potential is a linear function of vertex separation h 
for various quark- vertex distance L. To the right: It shows that penta-quark potential is a linear 
function of total flux tube length(solid line). The dashed line indicates a deviation from linearity 
in the UV region thanks to either the correction from attraction between vertices or the flip-flop 
mechanism, as also predicted in the lattice simulation. 

of N = 3, for instance, a baryon made of three quarks (QQQ) can combine with an anti- 
baryon made of three anti-quarks (QQQ) to form a tetra-quark (QQQQ), by annihilating 
a pair of (QQ) and reconnecting two baryon vertices. Similarly, a penta-quark (QQQQQ) 
can be constructed by an anti-baryon and two baryons. In general, one could form a bound 
state of (mN — 2m + 2) (anti-)quarks with m vertices as long as its energy is lower than 
sum of m (anti-)baryons 2 . From now on, we restrict ourselves to the case of iV = 3. It is 
straightforward to write down inter-quark potential of a tetra-quark, i.e. 

E 4Q = -4 j + (4L + h)a + 2C, (13) 

where h is the separation between two vertices. For small h, however, the lattice data tend to 
agree with the potential of two-meson, i.e. twice the potential given by ([6]) with separation 
h. This flux-tube recombination between connected four-quark state and two-meson state 
is referred to the flip-flop. The FIG. 2 shows that this flip-flop can also be realized in our 
model. Next, inter-quark potential of a penta-quark reads, 

E 5Q = -5y + (4L + 2h)a + 3C, (14) 

2 Here we simply mention those bound states made of single open chain of baryon vertices. In addition, 
there might be states with closed chain(s). We thank Kazuyuki Furuuchi to point out this. 
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FIG. 4: The solid curve shows the temperature dependence of string tension (normalized to unity). 
It can be very well fit by the dashed curve given by equation (|18p for low temperature. 

where we have assumed that three vertices are placed in equal separation h. From the 
equation above, it is obvious to expect linear relation between potential E and separation 
h as shown in the left plot of FIG. 3, or between E and total length of color flux tube, say 
4L + 2h, as shown in the right plot of FIG. 3. 



V. POTENTIAL AT FINITE TEMPERATURE 



Now we are ready to study the same potential at the finite temperature, which corresponds 
to introducing an AdS black hole, i.e. 

ds 2 = (^) 2 (-/(r)rft 2 + dx 2 ) + {-f^-dr 2 + R 2 dnl (15) 
it 'J v ) 

where f(r) = 1 — rf^/V 4 and the event horizon is related to the Hawking temperature by 
tt = irR 2 T. In the deconfmed phase where > R, the vertex falls into the black hole and 
leaves behind N dissociated quarks with total mass 3 , 

E = -N^^. (16) 
27ra' v ; 

However, the observers living in the boundary still see the vertex since it takes infinite time 
for the vertex to cross the horizon. In the confining phase where tt < R, one may either 
argue that black hole is unaccessible to us so the potential has no temperature dependence 



3 This bare mass is not physical but by convention of renormalization since ideal heavy quarks should have 
infinite mass which is unmcasurable. 



s 



for thermal AdS, or one may pretend that black hole still curves the bulk geometry. In the 



latter situation, the baryon potential has the following linear component 33| 




15_l 8p 4_ p 8 



(y 4 - p 4 )(16?/ 4 - 15 + 2p 4 + p 8 

where p = r^/r . Here, different from the situation at zero temperature, it is not allowed 
to send r — > 0. In fact, the reality condition requires that r > tt- In the FIG. 4 we 
see the temperature dependence of the effective string tension, which can be best fit in low 
temperature by 



(7(0) 



P 4 - (18) 



VI. IMPROVEMENTS 

We have argued in section HTT1 that our model can not completely fit the lattice result. In 
this section, we provide three possible improvements to the potential but also discuss their 
limitations. 



A. UV perfection 

As mentioned before that our model cannot completely fit the lattice result, introducing 
another UV cut-off r c may resolve this. However, this implies that instead of our original 
assumption of (infinitely) heavy quarks, the Cornell potential is realized via quarks with a 
finite mass 

If this is the case, we should also consider string breaking (pair production) as inter-quark 
separation grows. A careful examination excludes the possibility of sending ro — > 0, which 
leaves more than one choice for combination of ro and r c , or equivalently yn and y c = r c /ro. 
In fact, we only find reasonable solutions as yn < 2. One of these choices is that, for instant, 
yn ~ 1.50 and y c ~ 3.29. 
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B. IR refinement 



We start with a truncated model with IR cut-off at r = R. This implies that the 
confinement/deconfinement transition temperature is at T c = 1/tcR, which is ~ 212 MeV 
for the first fit A and ~ 312 MeV for the second fit B. However, this model may be refined 
with a different IR cut-off tir. This provides a new parameter for perfect fit, where 

C = ^r IR . (20) 

or 

with the choice of rm ~ 6.9891, However, this cut-off corresponds to the transition 
temperature [ljj , 

T c = ~ 54 MeV, (21) 

which is, unfortunately, far too low in comparison with the lattice prediction T c ~ 190 
MeV[38j. 



C. Correction via closed string channel and flip-flop mechanism 

The lattice simulation indicates that in the FIG. 3, relation between potential and length 
is not completely linear, especially for small separation h. The residual force outside each 
interacting group of quarks may help explaining the deviation from linearity. This could 
be realized in our model by two means: one way is to include the graviton (closed string) 
exchange between vertices. The other is the flip-flop mechanism in short h. Both give the 
expected 1/h behaviour in the IR region. 



VII. CONCLUSION 



In this paper, we have constructed static potential of heavy baryon and other exotic 
multi-quark configurations in the hard- wall model of AdS/QCD. In particular, we obtained 
a Cornell-like potential and make naive comparison with that in the quenched lattice 
calculation for SU(3) color groupf^. We found that the flux tube recombination between 
four-quark and two-meson states (flip-flop) happens when two vertex approaches each other 
in a tetra-quark configuration. In a penta-quark configuration, static potential is always 
linear with respect to inter-vertex distance or total length of flux tube, though a deviation 
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from linearity is expected for short h. We later proposed three different ways to improve 
our model. For UV perfection, we introduced a UV cut-off and there were more than one 
choices to achieve the goal. Nevertheless, the quark mass became finite and one has to 
take into account pair production (string breaking) as inter-quark distance grows. For IR 
refinement, we traded a tunable IR cut-off for perfect fit. However, this modification gave us 
far too low confinement /deconfinement transition temperature and was not so impressive. 
Finally, we argued that the deviation from linearity in the UV region in the figure fllllj) 
may come from correction via vertex interaction or flip-flop mechanism when the potential 
of two-meson state dominates. At end end, we would like to make a few remarks. First of 
all, there is no obvious reason to believe that AdS/QCD construction is trusted for large 
iV as well as in the AdS/CFT correspondence. On the top of that, we should not ask too 
much for quantitative accuracy even for small N, especially without taking in account the 
1/N correction. This may well explain that we did not find a perfect fit without causing 
another drawback. Secondly, in the confined phase, the string tension is exp ected to have 



a temperature correction cx — T 2 instead of equation (|18l) . observed by [l4j as well. It is 
interesting to see if we could settle this disagreement in the context of AdS/QCD. Thirdly, 
there was early debate on V-shape or A-shape flux tube configuration in the baryon ground 
state. Our construction always prefers the former since the latter costs more energy for 
slightly longer tube length and two extra vertex. At last, our proposal of UV perfection 
has similar effect when flavored quark is introduced 4 , and a warping geometry with linear 



dilaton is responsible to the string breaking 39] . It would be interesting to realize this 



breaking scenario in the soft-wall QCD model. In conclusion, the new feature as well 
as advantage of this construction, in comparison with the old version QCD string, is to 
realize all parts of Cornell potential, though not perfect, by one piece of configuration of 
strings and vertices stretching over truncated AdS bulk. This implies that there may be 
a chance in AdS/QCD construction for a unified theory of QCD in both perturbative and 
non-perturbative regions. 



The author would like to thank Carlos Nunez for pointing out this similarity. 
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